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Abstract 

We estimate the Tevatron Run II potential for top and bottom squark searches. 
We find an impressive reach in several of the possible discovery channels. We also 
study some new channels which may arise in non-conventional supersymmetry mod- 
els. In each case we rely on a detailed Monte Carlo simulation of the collider events 
and the CDF detector performance in Run I. 



1 Introduction 

For most of the next decade, the Fermilab Tevatron collider will define the high energy 
frontier of particle physics. The first stage of the Tevatron collider Run II, scheduled to 
begin in 2000, will deliver at least 2 fb _1 of integrated luminosity per experiment at 2.0 
TeV center-of-mass energy; this is more than 10 times the luminosity delivered in previous 
collider runs at 1.8 TeV. Major upgrades of the CDF and DO detectors are under way. 
Among other features, the detectors will have the ability to trigger on displaced vertices 
from bottom and charm decays using a precise microvertex detector. 

Along with top physics, searches for supersymmetry (SUSY) are among the main 
priorities for Run II. The Minimal Supersymmetric Standard Model (MSSM) is an at- 
tractive extension of the Standard Model. On the one hand, it agrees with precision 
measurements and cannot be easily ruled out through indirect searches [0. On the other 
hand, it is theoretically well motivated by various ideas about physics at very high energy 
scales - string theory, supergravity, grand unification. It also stabilizes the Higgs mass 
against radiative corrections, even though it does not provide the complete answer as to 
why the electroweak scale is so much lower than the Planck scale. 

Currently, we have no idea which of the supersymmetric particles are the lightest 
and as such, more easily accessible at colliders. With all things being equal, however, 
the lightest stop t\ is a very good candidate for studying at the Tevatron. First, it is 
colored, and so its production cross section is much larger than that for sleptons, for 
example. Second, it is feasible that the stop is the lightest squark. This may be due 
for example to a large mixing angle 9t between the superpartners of the left-handed and 
the right-handed top quarks, tt and Ir respectively, which form the lightest stop mass 
eigenstate: t\ = iicosOt + tRsin9 t . Alternatively, the large top Yukawa coupling X t 
enters the renormalization group (RG) equations of the stop soft masses and tends to 
reduce them in comparison to the other squarks. We shall discuss each of these effects 
in more detail in Section |2TT|. If the stop is the lightest of all squarks, its production at 
the Tevatron will be least suppressed by kinematics. Yet another motivation to look for 
a light stop is that it seems to be preferred for electroweak baryogenesis @. 

The purpose of this paper is to establish a basis for a systematic stop search in the 



upcoming Run II at the Tevatron. We shall consider the most promising stop signatures 
in a variety of supersymmetric models. In each case, we shall discuss under what cir- 
cumstances the stop can be light, what is the optimal search strategy, and what is the 
Tevatron reach. 

In general, the typical SUSY signatures are determined by the nature of the lightest 
supersymmetric particle (LSP) and whether R-parity is conserved or not. Conservation of 
R-parity implies that all SUSY decay chains end up in the LSP, which is stable and leaves 
the detector. If the LSP is charge- and color-neutral, as can often be the case in both 
supergravity and gauge-mediated models of supersymmetry breaking, then a typical SUSY 
signature is the missing transverse energy Ifix- A charged or colored stable LSP could lead 
to more exotic signatures (see, for example 0, [|, |J). R-parity violating interactions 
would allow the LSP to decay into Standard Model particles. For the purposes of this 
paper we limit ourselves to a class of SUSY models where R-parity is conserved and the 
LSP is charge and color neutral: i.e. it is either the lightest neutralino (as in supergravity 
mediated (SUGRA) models of SUSY breaking), or the gravitino (as in gauge mediated 
(GM) models). 

The plan of the paper is as follows. In Section [^ we start with a review of the traditional 
searches for a light stop, which were motivated to a large extent by the minimal version 
(mSUGRA) of the SUGRA models. Therefore, in Section |2~T| we first describe the relevant 
part of the mSUGRA parameter space which can be explored via those analyses. Notice, 
however, that even though we choose to work within the framework of mSUGRA, our 
results are valid for a generic MSSM. In fact, each analysis only makes an assumption 
about the mass ordering of a few supersymmetric particles, and is not constrained to any 
particular mechanism of SUSY breaking. In Section |2.2| we summarize the assumptions 
and basic facts from Run I that were used in the analyses to follow. We then proceed 
to estimate the Run II projections for the Tevatron reach in light stop searches for the 
following three channels: t — > xfb (Section |2.3| ), t — > \\ c (Section |2.4| ) and t — > i)£ + b 
or t — > £ + vb (Section |2.5|) . We then depart from the mSUGRA framework, and consider 
the possibility of non-universal scalar masses at the unification scale. This may lead to 
a situation where the neutralino LSP is mostly higgsino-like, which will alter the light 
stop search strategy. Such a case is studied in Section 0, where we present the reach as a 



function of the higgsino and stop mass. As the stop becomes increasingly degenerate with 
the LSP, the signals from direct stop production become lost and one has to look for stops 
among the decays of other sparticles, e.g. charginos. In Section ^ we consider decays of 
(possibly higgsino-like) charginos to stops. We devote Section § to sbottom searches, since 
in many respects the analysis is similar to some of the stop searches considered before. We 
delineate the relevant SUGRA parameter space and present Run II expectations. Finally, 
in Section || we translate our previous results for the case of light stops in GM models. 
We reserve Section |7| for our conclusions. 

2 Traditional Stop Searches 

In this section we discuss the relevant mSUGRA parameter space for the traditional stop 
searches. We then estimate the Run II sensitivity in these channels. 

2.1 Review of the relevant mSUGRA parameter space 

Most of the SUSY searches in the past have been performed within the SUGRA frame- 
work, where supersymmetry breaking, which takes place in a hidden sector, is commu- 
nicated to the MSSM fields through gravitational interactions. The sparticle spectrum 
can be calculated in terms of the soft SUSY breaking parameters, which are in princi- 
ple free inputs at the Planck scale Mp. In the minimal version of the SUGRA models 
(mSUGRA), it has become a custom to input the soft masses at the grand unification 
(GUT) scale Mqut instead, and in addition one assumes universality among the scalar 
masses, the gaugino masses and the trilinear Higgs-sfermion-sfermion couplings. The 
mSUGRA parameter space therefore consists of only five parameters: a universal scalar 
mass Mo, a universal gaugino mass M1/2, a universal trilinear scalar coupling term Aq, 
all defined at the GUT scale, as well as the ratio of the two Higgs vacuum expectation 
values v u /vd = tan/3 and the sign of the \x parameter. One can readily identify which 
regions of this mSUGRA parameter space would be associated with a light stop, simply 
by inspection of the stop mass matrix 

ka = ( M Qs + m 2 + {± -^ sin 2 6 W )M 2 Z cos 2(3 m t {A t + fi cot (3) \ 

* \ mt(A t + (i cot /3) M lR + m 2 + \sin 2 e w M 2 cos2(3 J ' [) 



where Mq (Mj ) is the soft mass parameter for the left-handed (right-handed) top squark, 
m t is the top quark mass, Mz is the Z-boson mass, 9w is the Weinberg angle and A t is 
the soft trilinear titjiH u coupling. All parameters entering eq. (|l|) are to be evaluated at 
the low-energy scale (near the stop mass). In the mSUGRA model, A t is most directly 
related to its boundary condition at the GUT scale, Aq, while Mq and M^ R have a 
dependence on both Mq and M1/2. The /1 parameter is obtained from the condition that 
proper electroweak symmetry breaking reproduces the experimentally observed Z-boson 
mass. Notice the possibility of a large stop mixing: 

• na 2m t {A t + ficot(3) 2m t A t ,.,.., , a \ fo \ 

smz&t ~ 5 ~ 7. — (at moderate to large tan p). (2) 

If there were no such mixing, the stop masses would be roughly equal to the corresponding 
soft mass parameters Mq 3 and M~ tR . The mixing, however, further reduces the smaller 
of the two mass eigenstates 0. (This effect is negligible for the first two generation 
squarks, because the corresponding quark mass entering eq. (H) is very small.) Since A t 
is directly related to its boundary condition Aq at the GUT scale, one would expect that 
mSUGRA models with large values of \Aq\ would have a light stop in their spectrum. 
As it turns out, a better parameter for light stop discussions is the dimensionless ratio 



a = Aq/ J Mq + 4M?/ 2 ~ A t I mi, which can be easily understood from eqs. (P and (|J). 
Increasing either Mq or Mu% would increase the diagonal entries in the stop mass matrix, 
and make the stop heavier, while Aq controls the size of A t at the weak scale. In Fig. [1] we 
show contours of the lightest stop mass (in GeV) in the (oo, tan/3) plane, for Mi/ 2 = 300 
GeV, Mq = 300 GeV and /i > 0. Inside the light-shaded region there is a scalar (stop, stau 
or CP-odd Higgs) which is too light or tachyonic. The dark-shaded region is theoretically 
allowed, but ruled out because of the CP-even Higgs mass limit of m^ > 95 GeV. In 
calculating the light stop mass, we have included the full one-loop corrections to the stop 
mass matrix || . This is needed in order to reduce the scale dependence of m^ , which is 
known to be significant, especially in such cases of large stop mixing. 

We see from Fig. [I] that a light stop is associated with a limited range of rather large 
values of |ao|, where the stop mixing becomes very large. We have checked, by scanning 
the mSUGRA parameter space, that if tan/3 > 2.0, values of ao between and 2 would 
never lead to stop masses below 200 GeV. 
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Figure 1 : Contours of the light stop mass rrif (in GeV) , versus the mixing parameter 
a and tan/3, for M 1/2 = 300 GeV, M = 300 GeV and fj, > 0. Inside the light- 
shaded region there is a scalar (stop, stau or CP-odd Higgs) which is too light or 
tachyonic. The dark-shaded region is theoretically allowed, but ruled out because 
of the CP-even Higgs mass limit of rrih > 95 GeV. 



Another possibility to have a light stop is that the soft mass parameters Mq and Mf 
entering eq. ([!]) are small by themselves. This may be due to the RG evolution down from 
very high scales, which tends to order the squark soft masses in an inverse hierarchy with 
respect to their Yukawa couplings. Since the top Yukawa coupling \ t is so large, the stops 
"feel" this RGE effect to a larger extent than the other squarks. As a result, the stop 
soft masses Mq 3 and Mi R are typically smaller than the soft mass parameters of the other 
squarks. This effect is strongest when the top Yukawa coupling is the largest, i.e. at small 
values of tan/3. Indeed, Fig. [I] exhibits a region at rather small values of tan/?, where 
one may find a light stop for a very wide range of the A parameter. Very low values of 
tan/3, however, are excluded both experimentally and theoretically. In the MSSM, the 
lightest Higgs mass m^ is correlated with tan/? and current bounds on m^ from LEP all 
but exclude values of tan/3 < 2.0. In Fig. we show contours of the light Higgs mass m h , 
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Figure 2: Contours of the light Higgs mass nih (in GeV) for the same parameters 
Fig. g[ 
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for the same parameters as in Fig. |T|. We see that the current Higgs bound already rules 
out the region where we get a light stop because of the large Yukawa RGE effects. This 
constraint may be relaxed in models with a non-minimal Higgs sector, like the NMSSM 
|J. In that case, there is still a theoretical lower bound on tan/3 < 1.3 — 1.5, which is due 
to the requirement that At remains perturbative when extrapolated up to very high scales, 
e.g. Mp. What is more, from Fig. |2| we also see that the direct Higgs searches in Run 
II |10[ also eat away from the light stop parameter space. If the Higgs escapes detection 
after the first stage of Run II (2 fb~ per detector), the expected rrih limit with (without) 
a neural network improvement of the analysis will be m^ > 120 (105) GeV |TD[. In that 
case, the mSUGRA light stop parameter space at a^ > — 1 (ao > 3) will also be ruled 
out. However, the light stop region at A < may actually be more effectively probed 
via stop searches rather than via the Higgs search. 

It is straightforward to check that the (very light) stop mass has no significant depen- 
dence on the other mSUGRA parameters. Our mSUGRA scan reveals that there exist 
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Figure 3: Contours of the light stop mass Mh (in GeV) in the Mq — Aq plane, for 
fixed Mi /2 = 300 GeV, tan/3 = 10 and fi > 0. 



very light stop solutions for almost any values of M and M\/ 2 \ The reason is that there 
always exist large enough values of \A \, which can increase the amount of stop mixing 
and yield an arbitrarily light stop in the spectrum. This is illustrated in Fig. [| where we 
show contours of the light stop mass in the M — A plane, for fixed M\i 2 = 300 GeV, 
tan/3 = 10 and \x > 0. This value of M\i2 is just beyond the Tevatron reach for chargino- 
neutralino searches in the trilepton channel |TI|| . We then see that for arbitrarily large 
Mq, i.e. heavy scalars, there is still a (rather limited) range of A , where of all SUSY 
searches in Run II, only the light stop searches have a chance of being successful. Of 
course, the branching ratios of the stop decays do depend on both M and My 2 , and this 
will determine the degree of applicability of each of the channels that we are considering 
below. We shall comment on the relevant range for M Q and M\i2 in each case as we go 
along. Also notice that the production cross section for stops is almost independent of 
any of the MSSM parameters, and is uniquely determined by the stop mass. 

We should point out that all stop search analyses below are pretty much model in- 



dependent and have only a few mild assumptions about the sparticle spectrum. The 
motivation for having a light stop in the spectrum, on the other hand, is very model- 
dependent. For example, in mSUGRA very large values of |a | violate the naturalness 
criterion fll2"f , implying that the odds for a light stop in mSUGRA are rather small. But 



the odds for mSUGRA being the correct model of SUSY breaking are probably pretty 
small too. So we adopt the stance that we should leave any theoretical prejudice behind 
and look for the stop in every conceivable channel, until the first experimental SUSY data 
are in. 

2.2 Experimental assumptions 

To determine the mass reach for various search modes in Run II we compare the number 
of expected signal events to a variation of the number of background events (assumed to 
be purely poissonian). A mass region is excluded if the signal is larger than 3 standard 
deviations of the background. 

The number of events is calculated as a product of the integrated luminosity, produc- 
tion cross section and efficiency after final cuts. 

The integrated luminosities of 2 fb _1 , 4 fb _1 and 20 fb _1 were used as expectations 
for different stages of Run II. 

At the Tevatron, third generation scalar quarks are expected to be produced in pairs 
via gluon-gluon fusion and quark-antiquark annihilation, and, consequently, the leading- 
order production cross section depends only on their masses. The next-to-leading order 
corrections increase the cross section and introduce a weak dependence on other masses 
and parameters (~ 1%) [[H|. Signal events were modeled using the PYTHIA generator 
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In the following we assume that the increase in the production cross section due to 
variation of the center-of-mass energy from 1.8 to 2.0 TeV will be approximately 40% for 
the squark pair production. Background cross sections were also scaled with factors deter- 
mined from the corresponding Monte Carlo generators [[0| (20% for the W/Z production 
and QCD). 

Our Run II efficiency estimates are based on the performance of the CDF detector in 
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Run I using the full detector simulation. When possible we extrapolated the results of 
existing CDF searches in corresponding channels. 



CDF is a general purpose detector described in detail elsewhere [I6|. The innermost 
part of CDF, a silicon vertex detector (SVX), allows a precise measurement of a track's 
impact parameter with respect to the primary vertex in the plane transverse to the beam 
direction. The momenta of the charged particles are measured in the central drift chamber 
which is located inside a 1.4 Tesla solenoidal magnet. Outside the drift chamber there 
is a calorimeter, which is organized into electromagnetic and hadronic components, with 
projective towers covering the pseudo-rapidity range \i]\ < 4.2. The muon system is 
located outside the calorimeter and covers the range in \rj\ < 1. 

Jet energies are calculated using the calorimeter energy deposition within a cone in 
f] — 4> space, where is the angle in the plane normal to the beam direction. The missing 
transverse energy is defined as the energy imbalance in the directions transverse to the 
beam direction using the energy deposited in calorimeter towers with \r]\ < 3.7. A lepton 
is identified by either hits in a muon chamber or a cluster of energy in the electromagnetic 
calorimeter, and an associated track in the central tracking chamber. 

Two heavy flavor tagging algorithms use the SVX information to tag charm and bot- 
tom quark jets, or c-jets and b-jets. In the Jet Probability ( JP) algorithm the probability 
that the track comes from the primary vertex is determined taking into account the im- 
pact parameter resolution. This probability is smaller for heavy flavor decay products 
because of their considerable lifetime. The track probabilities for tracks associated to a 



jet are combined into JP [PH . We associate tracks to a jet by requiring that the track 
is within a cone of 0.4 in rj — 4> space around the jet axis. Distribution of JP is flat by 
construction for light quark jets, originating from the primary vertex, and peaks at zero 
for heavy quarks. Since the JP is a continuous variable, the tagging can be optimized 
both for charm and bottom jets. In another heavy flavor tagging algorithm, SECVTX, 
a jet is identified as a bottom quark candidate if its decay point is displaced from the 



primary vertex [IB|. This algorithm is not very efficient for the charm tagging. 

In calculating Run II efficiencies we scaled by a factor of two the Run I efficiency 
for heavy flavor tagging. This factor accounts for the increase in the SVX geometrical 
acceptance due to the larger SVX length. Conservatively, we did not assume any other 
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improvements in the detector performance such as, for example, a better coverage for the 
lepton identification. The description of the upgraded CDF-II detector can be found in 

n. 



2.3 Reach in the bj£$T channel 

If the chargino (either gaugino-like or higgsino-like) is lighter than the light stop, the 



dominant decay of the stop is t — > bxt [P0j Plj p22[ til HI- If there are no sfermions 



(squarks or sleptons) lighter than the chargino, the latter decays to a real or virtual W 
and the lightest neutralino. In this case stop decays produce top-like signatures: 2Ws and 
2 b-jets. The only differences are kinematical : higher fa due to the massive neutralinos, 
different jet spectra and angular distributions. 

There are two possible search strategies for this decay mode based on different signa- 
tures : bijfa and £ + £~jfa. We note that in the case of the dilepton signature one pays 
a price of the low lepton branching ratio twice. Here we present a sensitivity study based 
on the bijfa signature. We select events with an isolated electron or muon with pt > 10 
GeV/c passing lepton identification cuts, and at least two jets, one with Et > 12 GeV and 
the second with Ej- > 8 GeV. At least one of the jets is b-tagged with the SECVTX al- 
gorithm. To decrease Drell-Yan and Z° background we removed events with two isolated, 
opposite sign leptons. For the reduction of QCD background we also required fa > 25 
GeV and A(p(fa - nearest jet) > 0.5. 

The main remaining backgrounds for this search are from the W + jets and top pro- 
duction. Table [l] lists the relative contributions of different backgrounds after all cuts 
for this channel together with two other experimental stop signatures to be considered 
later. In the table we quote the W(— > rz/)+jets contribution separately because of pos- 
sible hadronic decays of r which mimic a jet. The last line shows the total background 
cross section after final cuts. 

Typical values of efficiencies after final cuts for this and other (considered below) stop 
signatures are shown in Fig. |4]. The efficiencies are plotted as a function of the mass 
difference between the stop and its supersymmetric decay product. The mass difference 
effectively determines the kinematical properties of the reaction and, therefore, the ef- 
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Background 
process 


Stop signature 


U 3 $t 


cc$ T 


e+t-jfa 


W(-> e(/i)z/)+jets 


52% 


1% 


- 


W(-+ rz/)+jets 


3% 


52% 


- 


Drell-Yan,Z+jets 


2% 


8% 


26% 


wiy/wz/zz 


- 


3% 


13% 


« 


21% 


5% 


21% 


QCD (includes 66) 


20% 


23% 


41% 


Total cross section, fb 


980 


160 


50 



Table 1: Relative contribution of various backgrounds for the stop searches 
after final cuts and the total background cross sections. 

ficiency. Straight line fits show the efficiency parametrizations used for our estimates. 



We show the reach of Run II in this channel in Figure |5|. Even with 2 fb~ of data we 
will be sensitive to stop masses up to m t . The sensitivity vanishes when we approach the 
kinematic limit for this channel, because the b-jets become too soft. We see also that the 
projected reach for 2 fb~ completely overlaps with the region of expected sensitivity for 
chargino searches at LEP, but with higher luminosity the Tevatron will be able to extend 
those bounds. 

Note that the leptonic branching ratio of the chargino can be significantly increased, 
if the (left-handed) sleptons are relatively light, which happens when Mq < Mi/ 2 . This 
would lead to a much larger Tevatron reach in this channel ||25|| . In the mSUGRA model 
we find that for M < 0.55Mi/2 the electron and muon sneutrinos are lighter than the 
chargino, which would increase the chargino branching ratio to leptons to almost 100%. 



In this case the dilepton signature, discussed in detail in |2l| may give a better sensitivity. 

When the sneutrinos are light, the leptonic branching ratios of gauginos are high and 
this region of parameter space can also be effectively probed via searches for chargino/neutralino 



production in the clean trilepton channel [ll|. 
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Figure 4: Typical values of efficiencies after final cuts for the stop signatures con- 
sidered in the text. 



A possible improvement in the analysis would be a requirement of an additional b- 
jet, given a much higher integrated luminosity in Run II. This cannot be afforded for 
Run I analyses because the signal efficiency becomes too low despite the formally better 
sensitivity. 



2.4 Reach in the cctfr channel 

This is the simplest, and in some sense, most model-independent situation, which arises 
whenever the stop is the next-to-lightest supersymmetric particle. Then, the only two- 
body stop decay still open is t — > ex? 0, ^2], ^3] . In the absence of any flavor- changing 
effects in the squark sector (which is rather unlikely), this decay proceeds through a 
loop, otherwise it occurs at tree-level through stop-scharm mixing. The stop signature in 
this channel is two acolinear charm jets and missing transverse energy carried away by 
neutralinos. The events for this analysis in Run I were collected using a trigger which 
required J^t > 35 GeV. We select events with 2 or 3 jets with Et > 15 GeV and |r/| < 2. 
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Figure 5: Sensitivity of the light stop search in the Mjlfcr channel for several 
integrated luminosities. The analysis assumes 100% branching ratio for t — ► bxf 
and VFdike branching ratios for the chargino decays. 



The JpT cut is increased beyond the trigger threshold to 40 GeV and we require that the 
$t is neither parallel nor anti-parallel to any of the jets in the event in order to reduce 
the contribution from the processes where the missing energy comes from jet energy 
mismeasurement: min A(p(]fi T ,j) > 45°, Aip(]pT,ji) < 165°, and 45° < A(/?(j 1 ,j 2 ) < 
165°, where the jets are ordered in E T . We also veto electrons and muons to suppress the 
W+jets background. 

We use the JP algorithm to tag a charm jet requiring that at least one jet has a 
probability less than 0.05. This requirement, chosen to optimize the expected signal 
significance, rejects 97% of the background while its efficiency for the signal is 25%. 

As can be seen from Table [1|, the dominant source of remaining background for this 
analysis is VF+jets production where the vector boson gives a hadronically decaying r 
lepton. There is also a contribution from QCD multijet production. The middle plot 
in Figure |] shows the signal efficiency as a function of the stop and neutralino mass 
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tj -> c x? or t, -> b W x° 
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Figure 6: The same as Fig. [|, for the light stop search in the cc$t and bj£$T 
channels. The 

(i^bw+x°x). 



channels. The cc$t (bji^r) analysis assumes 100% branching ratio of t — ► c%§ 



difference. 

Figure |^ shows the Run II exclusion contours in this channel for several integrated 
luminosities. CDF Run lb results for this channel were presented in 



The reach in 



stop mass is determined by the accumulated statistics while the reach in neutralino mass 
depends on the efficiency of our selection cuts, where the most limiting is the $t cut, 
effectively fixed by the $t trigger threshold. This is what determines the gap between 
the kinematic limit and the excluded region. In Run II CDF will have a possibility to 



use the secondary vertex information at the trigger level [27|. Addition of a displaced 
track requirement may allow to lower the Ifi? trigger threshold to 25 GeV and, therefore, 
significantly extend the excluded region to the kinematic limit. Even with the high Ifi? 
cut, we see that already with 2 fb _1 the Tevatron will be able to probe regions well beyond 
the sensitivity of LEP searches. 

If the mass gap between the stop and neutralino masses is larger than m& + Mw, 
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the three-body decay t — > bW + x1 opens up. We use the signature bijlftr discussed in 
Section [O] to estimate the sensitivity for this kinematic region (see Fig. ^) . 

2.5 Reach in the jl + l~$T channel 

If the stop is lighter than the chargino, but heavier than any of the sleptons, then the 
three body decay modes t — > b£ + i)g and t — > b£ + U£ become dominant [^. [2^. |28fl . Such a 



situation may readily arise in the minimal SUGRA model. As we mentioned in Sec. |2.3| , 
for Mq < 0.55 • Mi/2, the sneutrinos are lighter than the chargino. In those cases, there 
almost always exist values for A$ (both negative and positive) which will bring the light 
stop mass in between rriy and m-+. 

Since in this case the stop leptonic branching ratio is high, we can afford to require 
two leptons in the final state. A 6-tag is not required in order to save on jet acceptance. 
We select events with two leptons P T {h) > 8 GeV/c and P T (h) > 5 GeV/c, E T > 30 
GeV and at least one jet with Et > 15 GeV. In order to suppress the bb background, we 
require the leptons to be isolated - the calorimeter energy sum in a cone of 0.4 in rj — <p 
space around both leptons should be less than 5 GeV. Further cuts on the angle between 
either of the leptons and the missing transverse energy reduce the background from jet 
mismeasurement: Aip(£, Ifir) > 20° and A<^(dilepton system, jfiz>) > 20°. The main back- 
grounds are Drell-Yan dilepton production, top, bb and QCD multijet production (see 
Table |1|) . The right plot in Figure f| shows the signal efficiency as a function of the stop 
and sneutrino mass difference. We show the reach of Run II in this channel in Figure [7]. 
Again, we see that the Tevatron will be able to go well beyond the LEP limits. 

At large values of tan/?, the same two effects that can lead to a light stop - the en- 
hanced RGE renormalization and the larger mixing, can also make the lightest tau slepton 
significantly lighter than the electron or muon sneutrino/slepton. (Note that because of 
the large stau mixing, now the chargino can couple to the lightest stau, so we are not 
limited to considering only the sneutrinos, which are typically heavier.) Therefore, there 
are larger regions in the mSUGRA parameter space where one encounters the hierarchy 
irif < m~+ instead. In the minimal SUGRA model this is possible for M < 5Mi/ 2 . In 



x 



that case, the three-body stop decay i — > br^v T may be dominant, giving rise to the 
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Figure 7: The same as Fig. ||, but for the light stop search in the b£ + £~ j$T channel, 
assuming that the stop decays as t — > bt^vg and £ — ► b£ + U£. 



signature bbr + T fa- Since the leptons from the tau decays are rather soft, it may be ad- 
vantageous to consider a signature where we replace one of the leptons with an identified 
tau jet: Mtj$ t P . 



3 Higgsino LSP 

So far we have been considering only minimal SUGRA models, which assume universality 
of the scalars and gauginos at the unification scale Mqut- Theoretically, however, this 
assumption is not very well motivated. If gravitational interactions are the mediator 
of supersymmetry breaking, then universality would naturally hold at the Planck scale 
M-p instead, and may get modified by whatever physics there is between Mqut and 
Mp. If there is grand unification, however, the GUT symmetry will preserve universality 
multiplet by multiplet. Thus one may expect that sparticles belonging to the same GUT 
representation still have identical soft masses at the GUT scale; while sparticles belonging 
to different GUT representations may have different soft masses. Within the framework 
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of a SUSY SU(5) GUT, this implies that the universal scalar mass parameter Mo is now 
being replaced by 4 scalar mass inputs at the GUT scale: M\o, the mass of the two up-type 
squarks, the left-handed down squark and the right-handed selectron; M 5 , the mass of the 
right-handed down-type squark and the left-handed slepton doublet; and M Hl (M H2 ) - 
the soft mass for the down- type (up- type) Higgs doublet. This type of model has become 
known as the non-universal SUGRA model. 

The non-universality in the boundary conditions for the scalar masses affects the SUSY 
mass spectrum in several ways. First, and most directly, it may change the ratios of various 
scalar masses, for example the left-handed and right-handed charged sleptons. However, 
the values of the scalar masses at low energies also depend on the gaugino masses through 
the RGE evolution, and the effects from any scalar mass non-universalities become diluted 
in the limit Mq <C Mi/ 2 . A true test of universality will therefore require measuring 
several squark or slepton masses to a very good precision, something which may only be 
accomplished at the next linear collider. 

There are also indirect implications of scalar mass non-universality The higgsino mass 
parameter //, which is determined from the condition of electroweak symmetry breaking, 
is sensitive to the soft mass spectrum of the Higgses and third generation sfermions. In 
the minimal SUGRA model, it turns out that typically the higgsino masses are quite a bit 
larger than the gaugino ones, and as a result, the LSP and the lightest chargino are mostly 
gaugino-like[]. In the non-universal SUGRA model, where the two Higgs soft masses are 
free inputs at the GUT scale, we often find regions of parameter space where \/j,\ < M\, 
and as a result, the two lightest neutralinos and the light chargino are almost degenerate 
and mostly higgsino-like. Relaxing the universality assumption for the gaugino masses 
may also lead to higgsino-like LSP PT| . 

The implicit assumption for all three stop searches in the previous Section was that 
the LSP is gaugino-like. In case of a higgsino-like LSP, the search strategy obviously needs 
to be modified. 



The analysis in Sections 2.3 relies on the presence of a hard lepton from the chargino 



1 There is one possible exception to this rule - for multi-TeV Mq one can find natural [[12) regions of 
|^| < M\,M% |3(J, where M\ and M2 are the SUSY breaking mass parameters for the U(l)y and SU(2)l 
gauginos, respectively. 



17 



tj -» b h{ or h{-> tj b -> X?(c)b 




150 200 



250 300 



M(tj),Gev/c 



Figure 8: Sensitivity to i — ► bh + or h + — ► bt decays for different integrated lumi- 
nosities. 



decay. If the chargino is higgsino-like, it is very close in mass with the LSP and the 
leptons from the chargino decays are too soft to be used for either triggering or off-line. 
The analysis in Section |2.5| assumes the existence of a significant mass gap between the 
LSP and the lightest chargino, bigger than the gap between the LSP and the stop itself. 
This is not true if the LSP is higgsino-like. In fact, the dominant stop decay in that case 
is t — > bxt • The subsequent chargino decay to Xi is associated with very soft leptons or 



jets. The only observable signature therefore is bb$T p2|. It is similar to the signature 



considered in Section 2.4 , except that now the heavy flavor jets are b-jets. To understand 



our sensitivity to stop in this channel we apply all the cuts discussed in Section |2.4j . We 
gain some sensitivity with respect to the cclpT analysis because the heavy flavor tagging 
technique that we use is more efficient to bottom than to charm. 

A family of curves below the diagonal in Fig. |8] represent the reach of Run II in the 
bbl^T channel for different integrated luminosities. We see that we restore our sensitivity 
to stop discovery even though the higgsino decay products are lost. The reach is similar 



to the one presented in Fig. |6] for the cclpT channel. 

4 Stop search in chargino (higgsino) production 

It may be that the stop is light, yet cannot be found by regular means since it is almost 
degenerate with the LSP and its decay products are too soft to be observed. Such re- 
gions of the SUGRA parameter space readily exist, and under those rather unfortunate 
circumstances one should look for stops among the decay products of some other, rather 
light particles^. Let us review our options. 

Of all the SM particles, the only ones likely to be heavier than the stop are the top and 
possibly the Higgs(es), but the production cross sections for the latter at the Tevatron are 
too small to be relevant. Top quarks may in principle decay to stops and gluinos [55]. This 



channel is usually closed, since the gluino is typically quite heavy. Some unconventional 
models [[| predict a light gluino in the tens of GeV range, with the gluino possibly being 
the LSP. However, existing data already rule out the range of gluino masses for which the 
two-body decay t — >■ tig is open ||. Top quarks may also decay to stops and neutralinos 
P2| . |5B| ]. One can look for these decays through precise measurements of the top branching 
ratios. If the stop is really degenerate with the LSP, it decays invisibly, and as a result 
the signature is an invisible top. If the t — Xi m ass difference is large enough so that the 
c-jets can be detected, yet small enough to evade the stop search in the cc$t channel, 
then the signature will be top quarks decaying to t — * c$t- The top cross section is big, 
but the width into Wb is also quite large, so these will be quite challenging analyses. 

We now turn to discuss the possibility of producing stops in SUSY cascades. Among 
the remaining SUSY particles, gluinos have the largest production cross section, and they 



can decay to tt pairs |21], |22], qTj. Since the stops are invisible, the signature is similar to 
the leptonic channels of top pair production. The crucial difference from tt production is 
that because of the Majorana nature of the gluino, half of the time the top quarks will 
have the same sign. Such an analysis is also in preparation for Run II. 

One can also consider neutralino decays to top-stop. In this case, however, the neu- 



2 In the case of an extreme degeneracy with the LSP (which is somewhat preferred on the basis of relic 
density arguments f33fl), the stops can be long lived and form bound states, "top squarkonia" |p4J] , which 
eventually decay through t\t\ annihilation. 
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tralinos would have to be heavier than 777.4 + m^. Since their cross sections are electroweak, 
they would be too small to be observed at the Tevatron. 

Stops may also appear in sbottom decays: b — > tW~ or b — > tH~, but first, these 
processes will have to compete with b — ► bx\ , which is preferred by phase space, and 
second, since stops are invisible, the final state signatures (VF+V4 7- , W + H~ or H + H~) 
will have very large backgrounds. 

Finally, we can consider production of charginos, which later decay to stops: x + ~~ * tb 
f32|j . This case looks more promising than the neutralino decays to stops. First, there is 
much less phase space suppression, and second, for gauginos, the chargino pair production 
cross sections are larger than the neutralino ones. In the rest of this Section we shall 
consider this channel in more detail. 

We start by assuming that the LSP is mostly Bino. (If it were higgsino, we are back 
to the case discussed in Section |3|). Then, if the SUGRA relations among the gaugino 
masses hold, the gaugino-like chargino would have to be twice as heavy. We shall consider 
fj, as a free parameter, which can be easily accounted for by non-universalities as discussed 
previously. Then, if |/i| <C M 2 , the lighter chargino will be higgsino-like; while if |/x| S> M 2 , 
we are back to the typical mSUGRA case of gaugino-like chargino. 

Since the stop is almost degenerate with the neutralino by our assumption in this 
section, the resultant charm jets from its decay are very soft and cannot be detected, so 
the final state signature will be bblpT and we can consider the same selection which was 
used in Section |2.4|. Yet the production cross section is different, since in this case we 



need to produce charginos via a weak process. In Fig. ^| we show a plot of the signal cross 
section times branching ratio versus the lightest chargino mass m^+ . The mass, as well as 
the gaugino/higgsino mixture of the lightest chargino in the figure are varied by changing 
/i, simultaneously adjusting A t so that to keep the lightest stop degenerate with the LSP. 
The rest of the sparticle spectrum is assumed to be very heavy. We have also included 
the contributions from relevant processes with heavier neutralinos which may decay to 
Xi- We use conservatively leading order results for chargino/neutralino production; next- 
to-leading order QCD corrections can increase the signal by 20-30% |38| . We see from 
Fig. |9] that the signal cross section is rather small, and for the range of chargino masses 
beyond the LEP coverage, the signal is only a few hundred fb. 
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Figure 9: a * Br = J2i,j a {XiXj )Br(xf —> hb)Br(Xj — > i*b) as a function of the 
chargino mass m-+, which is varied by changing \i. We fix Mi ~ 1M\ = 180 GeV, 
and appropriately adjust A t so that the lightest stop is degenerate with the LSP. 
The rest of the sparticle spectrum is taken to be very heavy. 



Our reach of Run II in this channel is presented by a family of curves above the 
diagonal in Fig. ^|. Again, we can recover some region in parameter space by using an 
alternative signature in our search, but the absolute reach is not very impressive, mostly 
because of the small production cross sections. 

5 Light bottom squarks 



Light sbottoms in the mSUGRA model can appear only at small Mi/ 2 and small M . In 
addition, they are always accompanied by light stops as well. In fact, throughout the 
whole mSUGRA parameter space, m~ h > m^. The only exception appears at values of 
tan/3 > 20 and \x < 0, where we find that —60 GeV < m~ b — mi < 0. (The correlation 
with the sign of /i is due to the SUSY threshold corrections to the bottom Yukawa coupling 
|55fl.) However, this part of parameter space is severely constrained j|DJ by the b — > sj 
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Figure 10: Sensitivity of the light sbottom search in the bottom neutralino channel 
for different integrated luminosities. 



measurement from LEP [[41]. Similar conclusions hold even when we relax scalar mass 
universality. 

In any case, having done the analysis for stops decaying to higgsinos, it is straightfor- 
ward to extend it to the case of sbottom production and sbottoms decaying directly to 
LSP's. Again we can consider the same selection used in Section |2.4j to estimate our sen- 
sitivity to a direct sbottom search. Our reach in Run II is presented in Fig. [H]. Although 
the cross section of sbottom production is roughly equal to the stop production cross 
section for the same squark mass, our reach in sbottom mass is somewhat higher than 
in stop mass in the similar cclpT channel - compare to Fig. |6|. This is due to the higher 
bottom tagging efficiency as compared to charm tagging efficiency. 

By scanning the mSUGRA parameter space, we found that the universality assumption 
leads to the relation m^o < 0.25 m~ h so that the mSUGRA parameter space maps onto 



the region below the dotted line in Fig. [10]. Relaxing scalar mass universality, we still 
find that typically m^o/mi < 0.25, but this ratio may go up to 0.4 for tan/5 > 20 and 
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/i < 0, which is in conflict with b — > 57. This means that the observation of a signal in 
this channel already in Run II, if interpreted as sbottom production, will hint towards a 
more unconventional low-energy SUSY, for example non-universal gaugino masses. 

6 Stop (sbottom) as NLSP in gauge-mediated mod- 
els 

Gauge mediation is an intriguing alternative for communicating supersymmetry breaking 
to the visible sector (MSSM) p2"| . It offers the potential of solving the supersymmetric 
flavor problem, and leads to novel collider phenomenology 



The minimal gauge-mediated models, where the only SUSY breaking contributions to 
the scalar masses are from SM gauge loops, do not predict a light stop in the spectrum. 

or 



However, one can easily imagine non-minimal extensions with extra gauge groups |l 1 



more complicated messenger sectors jDj), which may lead to a light stop. Gauge mediated 
models are characterized by a Goldstino LSP G, which is almost massless. Of course, all 
our previous results hold for the case of gauge mediated models with a stable neutralino 
NLSP, since in that case the phenomenology is no different from SUGRA models. But 
what is more, our results from Sections ^]4| and Section |5] may be applied for the case 
of stop and sbottom NLSP, which decays promptly to the Goldstino. Identifying the 
neutralino LSP with the Goldstino, and taking the limit m^o — > 0, we can read off the 



Run II stop (sbottom) mass reach from the x-axis in Fig. [^ (Fig. [T(|). For example, if 

we can 



the dominant prompt decay of the stop is the three-body mode t\ — > bW + G [[46] . 

see that already with 2 fb~ the Tevatron will be sensitive to light stop masses up to 160 

GeV. 

Prompt stop decays in gauge-mediated models can be expected only if the SUSY 
breaking scale is very low. Otherwise, the stops first hadronize in supersymmetric sbaryon 
or mesino states and then decay over macroscopic distances, leading to events with highly- 



ionizing tracks or displaced jets and large Ifir [0- It was recently pointed out, that in the 
case of prompt decays, stop mesino-antimesino oscillations can provide a very distinctive 
signature of like-sign top quark events p? . 
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7 Conclusions 

There are various ways to look for a light stop, depending on the rest of the sparticle 
spectrum. We presented a long list of stop signatures which can be looked for in Run II 
of the Tevatron with upgraded detectors. There are also numerous other possible stop 
signals. For example, models with broken R-parity allow for a set of stop decays beyond 
the ones considered here [J48| , as well as associated (single) stop production processes 
|49|. Also, the 4-body decay of the stop can become dominant in certain models with 



suppressed flavor-changing effects |Kj 



In most cases the upcoming Tevatron runs can extend the LEP-II reach [51]. 
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